Quantum mechanical cumulant dynamics near stable periodic orbits in phase space: Application to 
the classical-like dynamics of quantum accelerator modes 
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We formulate a general method for the study of semiclassical-like dynamics in stable regions of a mixed 
phase-space, in order to theoretically study the dynamics of quantum accelerator modes. In the simplest case, 
this involves determining solutions, which are stable when constrained to remain pure-state Gaussian wavepack- 
ets, and then propagating them using a cumulant-based formalism. Using this methodology, we study the rela- 
tive longevity, under different parameter regimes, of quantum accelerator modes. Within this attractively simple 
formalism, we are able to obtain good qualitative agreement with exact wavefunction dynamics. 
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I. INTRODUCTION 

Since their initial discovery quantum accelerator modes 
have proved to be a fascinating example 
of a robust quantum resonance effect, and an exciting devel- 
opment^ a^-optical studies im H 13 H H nil H [13 
[nl[^, Ts", 14] of quantum-nonlinear phenomena fl^]. The 
demonstrated coherence of their formation [4} promises im- 
portant applications in coherent atom optics 13 [3 [Si- The 
experimental configuration in which they have been observed 
is closely equivalent to those used by t he g roup of Raizen 
10, [13 lllll - and subsequently by others I12I in the study of 
quantum 5-kicked rotor dynamics, in particular in observing 
the prediction of dynamical localization ] 17]. In a configura- 
tion consisting of a laser cooled cloud of freely falling cesium 
atoms subjected to periodic 5-like kicks from a vertically ori- 
ented off-resonant laser standing wave |[ll0|10,|l|3l3lil, 
quantum accelerator modes are characterized experimentally 
by a momentum transfer, linear with kick number, to a sub- 
stantial fraction (up to ~ 20%) of the initial cloud of atoms. 
The dynamical system, including the explicit presence of 
gravity, we term the 5-kicked accelerator ]3]. The experimen- 
tal observation of quantum accelerator modes lead to the pio- 
neering formulation by Fishman, Guarneri and Rebuzzini ] 18] 
of the strongly quantum-mechanical dynamics in terms of an 
effective classical map. This is justified by the closeness of 
the kicking periodicity to particular resonant times. The re- 
sulting limiting dynamics are termed e-classical llT8l[l9il . and 
have been used to great effect in the interpretation and predic- 
tion of ejcperimentally observable quantum accelerator modes 

eiili- 

The 5-kicked accelerator is therefore also attractive in that 
it is possible to tune its effective classicality in an accessi- 



ble regime far from the true semiclassical limit, making it an 
ideal testing ground for semiclassical theories. Semiclassical 
approaches in quantum chaotic dynamics have proved very 
successful in forging conceptual links between classically 
chaotic systems and their quantum mechanical counterparts 
(T?]. When trying to include quantum mechanical effects, an 
obvious step beyond point-particle dynamics is to consider the 
evolution of Gaussian wavepackets. Straightforward semi- 
classical Gaussian wavepacket dynamics are limited in that, 
e.g., the wavepacket is unrealistically forced to maintain its 
Gaussian form. Pioneering work by Huber, Heller, and Lit- 
tlejohn ]20] proposed remedying this by allowing complex 
classical trajectories. These also permit the study of a wider 
range of classically forbidden processes, and the propagation 
of superpositions of Gaussians. We propose an alternative ap- 
proach, which, most simply, is to follow the dynamics of the 
cumulants of initially Gaussian wavepackets. When taken to 
second order, the dynamics are described purely in terms of 
means and variances, as in a Gaussian wavepacket, but evolu- 
tion into non-Gaussian wavepackets is not prevented. In this 
paper we develop an appropriate general formalism, and ap- 
ply it to the phenomenon of quantum accelerator modes, thus 
reaching beyond an e-classical description of the dynamics. 

The Paper is organized as follows: Section II details in a 
general way the necessary essential formalism on Gaussian 
wavepackets and non-commutative cumulant hierarchies, as 
well as specifying important differences between Gaussian 
wavepacket dynamics and second-order truncations of the cu- 
mulant hierarchy. Section III introduces quantum accelerator 
modes, as appearing in an atom-optical i5-kicked accelerator. 
This is followed by a comprehensive derivation of the kick- 
to-kick operator dynamics, crucial for the second-order treat- 
ments that follow. The section culminates with the derivation 
of the relevant e-classical dynamics, along the lines of Fish- 
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man, Guarneri, and Rebuzzini fTs*!. Section IV builds on the 
operator dynamics derived in the previous section to produce 
an approximate second-order cumulant description, in which 
the e-classical theory is a first-order expansion within the cu- 
mulant hierarchy. There follows a worked example of the nec- 
essary methodology to determine approximate periodic orbits 
within the resultant second-order kick-to-kick map, consider- 
ing the most significant quantum accelerator modes 1 1 ] . These 
approximate solutions are propagated with the second-order 
mapping equations and with the exact time-evolution operator, 
which yields useful insight into the experimentally observed 
finite lifetimes of quantum accelerator modes, demonstrably 
showing the utility of our second-order approach. Section V 
consists of the conclusions, which are followed by four tech- 
nical appendices, which serve to make the Paper entirely self- 
contained. 



II. OVERVIEW OF THE ESSENTIAL GENERAL 
FORMALISM 

A. Gaussian wavepackets 

We consider two conjugate self-adjoint operators: | and ^, 
such that [|, I] = irj, and a Hamiltonian described in terms 
of these operators H{^, I). The dynamics of these operators 
can be fully described by the expectation values ju^ = (|) and 

- {0 only as ?/ — » 0. In this limit there is a well-defined ^, 
^ phase space, which generally consists of a mixture of stable 
islands based around stable periodic orbits, and a chaotic sea. 
This is the case for the specific system we consider, the 6- 
kicked accelerator (see Fig.flJ fs']. 

When considering dynamics near a stable periodic orbit 
in phase space, we note the facts that: local dynamics ap- 
proximate those of a harmonic oscillator f2l'l, and Gaussian 
wavepackets remain Gaussian when experiencing harmonic 
dynamics. This can be used as a motivation for the initial 
use of a Gaussian ansatz of the form 1 20. .221 



X exp 



[1 - /2cr|^/?7][^ - /^J^ _ ^ 



4cr| 



(1) 



where cr| = (l^) - {^}^ is the variance in |, and cr|^ - {^l + 

(i)/2 - is the symmetrized covariance in | and ^. As 

Eq. Q describes a minimum uncertainty wavepacket, the ^ 
variance, crl - {^-} - {0-, can be deduced from the general 
uncertainty relation 



2 2 / 2 \2 



(2) 



This can be seen from Eq. Q, using -irjd/d^ as the ^ repre- 
sentation of If the stable islands around the periodic orbits 
of interest are significant compared to the size of a minumum 
uncertainty wavepacket, we may find stable periodic orbits in 
fi^, fj.(, cr|, and cr|^) when such a Gaussian ansatz is enforced. 
In general this stability can only be approximate, but we will 



nevertheless utilize such solutions, as they are good estimates 
to maximally stable Gaussian wavepackets. We note that it 
is possible to significantly extend semiclassical techniques to 
obtain good error estimates for Guassian evolutions, includ- 
ing for classically chaotic situations 123,1 . Such calculations 
require substantial sophistication, and our intent here is some- 
what different, in that we wish to determine useful information 
within a very simple description. 



B. Non-commutative cumulants 

A complete picture of the observable dynamics can only be 
determined from the time-evolution of all possible expecta- 
tion values of products of the dynamical variables. Except for 
very simple systems, this produces a complicated hierarchy of 
coupled equations. 

In order to gain any insight we must determine a truncation 
scheme to reduce this to a managable description. This is in a 
sense achieved by the Gaussian ansatz, which considers only 
means and variances. Means and variances are the first two 
orders of an infinite hierarchy of cumulants |24], which we 
denote by double angle brackets to distinguish them from ex- 
pectation values. The non-commutative cumulants can be ob- 
tained directly in terms of operator expectation values through 

m 



{{qi ■ --qn)) 



OTl dT„ 



(3) 



Ti=0,...,T„=0 



where e {|, 0. More conveniently, the expectation values 
can be expressed in terms of cumulants: 

(qi) 

iqiqi) =«§i§2» + «§i»«§2», 
{qmh) ^{{qmh)) + «§i»«§2§3)) + «§2))«§i§3» 
+ mrnqm)) + mmqimh)). 



where the ordered observables have been partitioned in all 
possible ways into products of cumulants. 

Cumulants tend to become smaller with increasing order, 
unlike expectation values or moments, which always increase. 
Intuitively, higher-order cumulants encode only an "extra bit" 
of information that lower-order cumulants have not yet pro- 
vided. It is therefore often possible to provide a good de- 
scription by systematically truncating, expressing moments 
of all orders in terms of cumulants up to some finite order 
j^S*, '26]. Truncating at first order is equivalent to consider- 
ing only mean values, and thus reproduces the corresponding 
Hamilton's equations of motion. 

It is tempting to think that truncating at second order is 
equivalent to enforcing a Gaussian ansatz. As will be shown 
explicitly in our model system, the 5-kicked accelerator, this 
does not in general reproduce the dynamics given by enforc- 
ing a Gaussian ansatz. Gaussian wavepacket dynamics are 
unitary, meaning that the uncertainty relation of Eq. (0 is al- 
ways exactly observed, and that one need consider only two 
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of (t|, cr^, and cr|^. This is only true when no terms in the 

Hamihonian are of greater than quadratic order in | and I 
j20]. Furthermore, finding a fixed point of /i^, yU^, cr^, and 

cr^^ is equivalent to finding a perfectly Gaussian eigenstate of 
the system, which can only be true for the harmonic oscillator. 

When propagating the second-order truncated equations of 
motion for the first and second order cumulants, it is necessary 
to consider the dynamics of each of cr|, cr^, and cr|^ explicitly, 
as the uncertainty relation is not explicitly hard-wired into the 
formalism. This implies that the evolution described solely in 
terms of the first and second order cumulants is not unitary. 
This feature of this approach more accurately reflects the fact 
that truncating generally leaves us with an incomplete descrip- 
tion of the dynamics, with a correspondingly inevitable loss of 
information about the state of the system. We are not, in prin- 
ciple, restricted to initially pure states, although this flexibility 
is not exploited in this paper 

Nonetheless, when situated inside a stable island in ^, f 
phase space, such a "stable" Gaussian wavepacket should be 
long-lived due to the harmonic nature of the local dynamics 
12711 . We can then use the equations of motion appropriate to 
second-order cumulant dynamics to get an idea of how long- 
lived the initial wavepacket actually is, as physically sensible 
imperfections are included in the dynamics in a straightfor- 
ward manner. 

The approach we have outlined is most obviously applica- 
ble in the standard semiclassical regime, but is not restricted 
to it. We will illustrate our method by applying it to a very in- 
teresting and experimentally relevant system, the quantum 5- 
kicked accelerator |3], outside the semiclassical regime (see, 
also, AppendixlAl. Our approach provides useful insights on 
the longevity of quantum accelerator modes in this system, 
essential for their possible application in coherent atom optics 



III. QUANTUM ACCELERATOR MODES 
A. Introduction to the (5-kicked accelerator 

1. Experiment 

In the Oxford experimental realization of the quantum 5- 
kicked accelerator, ~ 10^ cesium atoms are trapped and 
cooled in a MOT (magneto-optical trap) to a temperature of 
5juK, yielding a Gaussian momentum distribution with a full 
width half maximum of 12 photon recoils j3- The atoms 
are then released and exposed to a sequence of equally spaced 
pulses from a standing wave of higher intensity light 15 GHz 
red-detuned from flie 6^5 1/2 ^ 6^^1/2, (f = 4 ^ = 3) Dl 
transition. 

In this parameter regime, the spatial period of the stand- 
ing wave is /Ispat = 447 nm, and the half-Talbot time Ti /2 - 
66.7yus. The half-Talbot time is equal to the pulse periodicity 
at which a quantum antiresonance would be observed if the 
initial condition were a zero-momentum plane wave |3, 
and is so-named in analogy to the Talbot length in classical 



optics f?, "2^. This quantity is of central importance, as it 
is when the pulse periodicity approaches integer multiples of 
T\i2 that quantum accelerator modes are observed experimen- 
tally 01 B II S II nil- The peak intensity in the stand- 
ing wave is ^ 5 X lO'^mW/cm^, and the pulse duration is 
tp - 500 ns. This is sufficiently short that the atoms are in 
the Raman-Nath regime for the observed momentum scales, 
and hence each pulse is a good approximation to a (^-function 
kick. The potential depth is quantified by 0^ - Dr'tp/^di, 
where Q. is the Rabi frequency and 6l the detuning from the 
Dl transition. 

In the Oxford experiment, the widths of both the beam and 
the initial laser-cooled atomic cloud are ~ 1 mm, and situa- 
tions exist where the consequent inhomogeneity of the driv- 
ing strength must be taken explicitly into account |4]. It is 
frequently sufficient to consider an averaged value only, how- 
ever Ol, which is the approach we adopt in this Paper This 
is also appropriate for a configuration where the initial atomic 
cloud is much smaller than the laser beam profile. 

After the pulsing sequence, the atoms fall through a sheet 
of laser light resonant with the 6^S 1/2 — > 6^^3/2, (F = 4 — > 
F" = 5) D2 transition, 0.5 m below the MOT. By monitoring 
the absorption, the atoms' momentum distribution can then 
be measured by a time-of-flight method, with a resolution of 
2 photon recoils. For more complete details of the Oxford 
experimental configuration, see Refs. |2l|3(]. 



Model Hamihonian 



The dynamics of the atoms in the Oxforc 
or modes experiment 



; Oxford quantum acceler- 
ator modes experiment |r,'2V3V4V':?,'^,'7'l are well modelled 
by the one-dimensional 5-kicked accelerator Hamihonian: 



^ = £^ + - + cos(Gz)] 2 5(f - «r) . (5) 



Here z is the position, p the momentum, m the particle mass, g 
the gravitational acceleration, t the time, T denotes the pulse 
period, G - In/Aspat where /Ispat is the spatial period of the 
potential applied to the atoms, and fi^^ quantifies the depth of 
this potential. 

As this Hamihonian is periodic in time, the time evolution 
of the system can be described by repeated application of the 
Floquet operator 



F - exp — 



— +mgz 
Zm 



r exp(j0rf[l -Hcos(G£)]), (6) 



which describes the time evolution from immediately before 
the application of one kick to immediately before the applica- 
tion of the next. 



3. Resonance condition 

The near-fulfilment of the quantum resonance condition 
(closeness to particular resonant pulse periodicities 
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|3Clj) means the free evolution of a wavefunction, e.g. initially 
well localized in momentum and (periodic) position space im- 
mediately after it experiences a kick, causes it to rephase to 
close to its initial condition just before each subsequent kick. 

The innovative treatment due to Fishman, Guarneri, and 
Rebuzzini accounts for this in terms of a so-called e-classical 
limit iTs", [l^ I30I1 . where a kind of kick-to-kick classical 
point dynamics is regained in the limit of the pulse period- 
icity approaching integer multiples of the half-Talbot time 
T112 = 2mn/nG^ \2], i.e., as e = 2n{T IT\i2 - ^) — > 0, where 
^ G Z. This accurately accounts for the observed acceleration 
for up to ~ 100 kicks, as well as predicting numerous ex- 
perimentally observed high-order accelerator modes |5]. It is 
this e, whose smallness indicates nearness to special resonant 
kicking frequencies, leading to the production of quantum ac- 
celerator modes 1 1 8 ] , and nol ti (or equivalently the commonly 
used dimensionless rescaled Planck constant k fioll . as defined 
in Appendix IaI. which takes the place of 77 in our cumulant- 
based approach. 

We now describe the treatment of Ref. ifTsIl to justify the ap- 
propriate phase-space which is the starting point of our anal- 
ysis, in a somewhat different, more operator-oriented form, 
which turns out to be more convenient for our purposes. We 
provide enough detail for the explanation to be self-contained, 
as well as setting notation and a conceptual basis for our 
second-order cumulant treatment, as applied to quantum ac- 
celerator mode dynamics. 



B. Derivation of the Heisenberg map 

1. Gauge transformation 

Moving to a frame comoving with the gravitational acceler- 
ation [Uit) - expiimgzt/K)], the Hamiltonian of Eq. (|5} trans- 
forms to U(t)HU^(t) - mgt We write the resulting gauge- 
transformed Hamiltonian as 



The time-evolution operator from f = to f = «r is then 
given by 



(9) 



T^n - F„F„-i . ..Fx -T 
where T" is a time-ordering operator. 

2. Block theory 



In this gauge the Hamiltonian is spatially periodic, and we 
can therefore invoke Bloch theory llSllSlll . Specifically, we 
parametrize the momentum eigenstates as \{hGY^p = k+ fi), 
where k e Z, fi e [0, 1), and /3 is termed the quasimomentum. 
The Hamiltonian's spatial periodicity implies that only mo- 
mentum eigenstates differing by integer multiple of fiG will 
couple together, i.e., the quasimomentum is conserved. 

Correspondingly, the position eigenstates are parametrized 
as \Gz - 2jil + 0), where I e Z, e [0, 27r). The position and 
momentum operators can be similarly parametrized as 



Gz = 27rl + 6, 



defined such that 



l\Gz = 2nl + 6}^ l\Gz = 2;rZ + 6), 
0\Gz = 2nl + 6)^ 0\Gz = 2nl + 6), 

k+P)^ mtiGT^p = k+P), 



k\(nGr^p 



mGT'p = k+/3}^ mGT'p ^k+p). 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 



Substituting Eqs. dlOt and ( II U into Eq. the Hamiltonian 
in the frame accelerating with gravity can therefore be written 



H = 



[nG{k + p)-mgtf . V T^ nA^ 
— n^d[^+cos{e)] d(t-nT), (16) 

H=-CO 



H 



(p - mgty 

2m 



n(Pd[l + cos(Gz)] Yj - «^)- (7) 



Due to the explicit time-dependence in the free-evolution 
part of the Hamiltonian, it is no longer periodic in time. 
The Floquet operator of Eq. which transforms to 

U{nT)FU'\{n - l]r), consequently also has an explicit time 
dependence, and it becomes necessary to specify which kick, 
and subsequent free evolution, it is describing. 

We therefore write the gauge-transformed kick-to-kick 
time-evolution operator, from t - {n - \)T io t - nT, as 



F„ =exp(-- 



—T-—(2n-l)T^ 
2m 2 



(8) 



X exp (i(pd [1 -I- cos(Gz)]) , 



where we have neglected an irrelevant global phase equal to 
exp(-/mg2[3„2 _ 3„ + l]T^/6n). 



where from the commutation relations (see Appendix |E] for 
derivations) [0, 6] - 0,k] - 0, we deduce that 0, H] - 0, and 
thus confirm that the quasimomentum is a conserved quantity 
in the accelerating frame. 



3. Time-evolution operator 



Substituting Eqs. JlOt and il l\ into Eq. (|8}, we get 



Fn =exp -/ 



(k^ + 2kp)- 



TiG^T 
2m 



- k{2n - 1) 



gGT^ 



X exp -i 



,nG^T 



2m 

X exp[i^j [1 + cos(0)j) 



-P{2n-1) 



gGT^ 



(17) 



where the middle, quasimomentum-dependent term com- 
mutes with the the other terms of the kick-to-kick time- 
evolution operator, and with all subsequently applied such 
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operators. After n kicks, the result is a phase deter- 
mined entirely by the initial quasimomentum equal to 
ex^{-in{0^nG^T 12m - pngGT^ 12]). 

Removing this relatively trivial time evolution from explicit 
consideration, recall that we are interested in the case whre T 
is close to integer multiples of the half-Talbot time Ti/a - 
2nmlhG^. Thus, substituting m T - {2nl + e)m/ftG^, and 
using that e 



-ik-nt 



F„ - exp \-i\—e + 2nk 



{ . T / 1 , 



(18) 



:exp(/0<i[l +cos(0)]), 



where we have inserted 27rf2 - gGT^, a dimensionless param- 
eter accounting for the effect of the gravitational acceleration, 
also known as the unperturbed winding number fs']. 

The total time-evolution operator from f = to f = nT is 
then given by 



exp -imi 



exp^-iTrn 

-pnQ. 

I 111 



1/2 



T 



(19) 



4. Effective Hamiltonian 

The dynamics governed by the time-evolution operator of 
Eq. (I18> are, just prior to each kick, exactly equivalent to those 
determined by the effective Hamiltonian 



the Heisenberg equations of motion for the dynamical vari- 
ables can be readily determined from the effective Hamilto- 
nian given in Eq. ( I20> . and solved to give a discrete kick-to- 
kick Heisenberg map: 



k„+\ = k„ - 4>d sin(0„), 
= hi + {k„^\ - 



- -|n + ^|Q+j6„ — 
2 \ 2) 



--\n+ -\q. + B„ — 

2 \ 2) ^ Ti/2\ 



(25a) 
(25b) 
(25c) 
(25d) 



Using this procedure to determine Eq. i25i is exactly equiv- 
alent to determining the kick-to-kick time-evolution directly 
from the time-evolution operator of Eq. ilHi by 0„+i = 
n+i^"^«+i- Equation (|533) is simply a confirmation that the 
quasimomentum is a conserved quantity, and from Eqs. ( 125 al l 
and (I25b> we see that the evolution of the conjugate pair 6 and 
k is completely decoupled from that of the discrete position 
variable I. Both of these properties are a direct consequence 
of the spatial periodicity of the gauge-transformed Hamilto- 
nian, combined with the (typical) chosen parametrization of z 
and p into discrete and continuous components. 



H = ®("^ - t)H„&{t -[n-l]T) 

oo 

- Hcpdll + cos(e)] ^ 6(t - nT) 



(20) 



Partial trace 



where is a step function. Essentially this means that the free 
evolution Hamiltonian 



^" " r 1 2"^ ^ 



1/2 



(21) 



can be considered time-independent between two successive 
kicks, but from before to after application of each kick, there 
is a stepwise change to the term linear in k. 



Corresponding to what in Ref. |18] is termed a Bloch- 
Wannier fib ration, it is convenient to exploit the fact that sub- 
spaces associated with distinct values of the quasimomentum 
P are decoupled, by considering the evolution of the relevant 
dynamical variables 9 and k, conditional to a specific value of 
the quasimomentum p, separately. 

Starting from a general density operator 



5. Heisenberg map 

Using the commutation relations (see Appendix |B] for a 
derivation) 



[0, /t] = i, 

[k,l] = [Vp] = ^, 
2ti 

\B,p\^m\^\Kp\^^, 



(22) 
(23) 
(24) 



: rr dpdp' dkk'(j3,p') 

X \inGr^p ^k+ P){{nG)-^p ^k' + p'\ 



(26) 



we consider expectation values determined by a form of par- 
tial trace, over the density operator multiplied by the dynami- 
cal observable. In general, we define 



6 



J CO 



(27) 



r 



where NiP) - Tj'k'=-oodkk(J3,fS) is a normalizing constant, 
causing the expectation value to be independent of the pro- 
portion of the population occupying any given /5 subspace. 
If we define the projection operator 



P(j3)^ ^ \(nG)-'p^k- 



k+Pl (28) 



which is that component of the identity specific to a particular 
quasimomentum subspace, then we can define the expectation 
values of interest by taking the full trace over the density op- 
erator multiplied by the operator of interest and a projector; 



we remove the explicit time-dependence and jS-dependence of 
the free evolution, to arrive at 

0„+i (y8) ^9„(J3) + sgn(e)J„^i(J3), (34a) 
J„MI3) =J„(/3) - ~k smCe„)P(/3) - ^gn{e)2nQ.P{fi) . (34b) 

With an appropriate definition of of the inner product, the 
dynamics can be mapped onto those for a rotor, the so-called 
y6-rotors of Ref. 1 18] (see Appendix|DJ. 

C. e-classical dynamics 



where Tx{pr{fi)} = N(J3). 



(29) 



7. Projected Heisenberg map 

Equation j29t shows that expectation values, which are in 
general time-dependent, of the operator-valued observables 6 
and k, multiplied by the projection operator P(j3), are exactly 
equal to the values determined by taking the y6-dependent par- 
tial traces of 6 and k, considered on their own (i.e., without a 
projection operator). 

It is therefore clear that the dynamics of and k, conditional 
to a particular value of the quasimomentum /3, can be equiv- 
alently described by the evolution of the projected operators 
0(fi) = ffPifi) and k(fi) = kP(/5). 

We can thus speak of a projected Heisenberg map: 



+ In 



{ 1\ T 



I„+i(/3) ^I„(J3) - ksm(0„)'P(J3), 

where we have introduced J = \e\k and k ■ 
ing implies the following commutators: 



(30a) 



(30b) 



\e\4>d- This rescal- 



[e,i] =/ki, 



(31) 
(32) 



which will therefore vanish as e — > 0, i.e., as T approaches 
rational multiples of the half-Talbot time Ti/2- 
If, furthermore, we introduce 



Jn(J3) = In(/3) + sgn(e)2;r 



1/2 



(33) 



1. Derivation oftiie e-classical map 

To determine the dynamics of the mean values (conditioned 
to a particular value of the quasimomentum fi), one simply 
takes the appropriately normalized expectation values of both 
sides of Eq. ( I34> . 

Due to the presence of sin(0„)!P(y6), this generally couples 
to an infinite hierarchy of higher-order moments or cumulants. 
However, in the limit e — > 0, we can make the substitution 



Tr{p0"!P(y3)) 
Tr{pf'(/3)) 



(35) 



Essentially the limit e — > allows us to treat the 6{fi), j'(J3) 
dynamics as those of a point particle, and higher order cor- 
relations can be discarded. This truncation is discussed more 
comprehensively in Section lTV Al 

The result is that the dynamics reduce to the kick-to-kick 
evolution of a pair of coupled scalar quantities: 

On^m ^0,m + sgn(e)J-„+i(;6), (36a) 
Jn^i (J3) =J,m - k sin 0„{li) - sgn(e)27rQ, (36b) 

where Oifi) and J'ifi) can be considered a pair of canonically 
conjugate classical action-angle variables. 

Quantum accelerator modes are explained by stable peri- 
odic orbits in the 0(j6) J^ijJ) phase space 1 18|. It is important 
to note that the dynamics described by Eq. j36> are indepen- 
dent of the value of y6, and that this is the reason we only need 
consider one phase space. We retain the fi argument due to the 
fact that if one begins with an initial cloud of atoms where all 
P subspaces are populated, then in general the inital conditions 
propagated by Eq. ( I36> are different, in a yS-dependent way. It 
is therefore necessary to distinguish these parallel evolutions 
if one wishes to finally produce a physically meaningful an- 
swer 

The true dependence of the dynamics on the value of /? has 
been removed by Eq. (I33> . which changes the frame in which 



the dynamics are observed from one falling freely with grav- 
ity, to one which is explicitly dependent on the quasimomen- 
tum p. To determine the dynamics of a cloud of falling atoms, 
it is therefore necessary tranform back the results of Eq. ( 13 6> 
using Eq. ( I33> . so that all quasimomentum subspaces are once 
again on an equal footing, i.e., simply falling freely with grav- 
ity. 

Measurements in the Oxford experiment are thus more 
closely related to the I{j3) distributions, the evolutions of 
which are explicitly y6-dependent [see Eq. (I30H . than the J"(y6) 
distributions, as will be explained in more detail in Section 
imC3lll8tl . 



2. Periodic orbits in the e-classical phase space 

The Oip), J^iP) phase space structure described by Eq. ( I36> 
is 2n periodic in J^ifi), meaning that although the phase space 
is cylindrical [cyclic in the Oifi) angle variable and infinite in 
the J^iP) action variable], it can be divided up into structurally 
equivalent phase-space cells. In long in the J^ifi) direction. 

We now consider an initial condition - J', 6(J3) - 

9, that lies exactly on a periodic orbit in phase-space. By 
definition the momentum of such an initial condition evolves 
as 



fa1 



Jnp(P)^J + 27Tnl 



(37) 



where p is the order of the periodic orbit (the number of 
points belonging to the periodic orbit contained within a sin- 
gle phase-space cell), and j is the jumping index [the number 
of phase-space cells traversed after p iterations of Eq. (I36H 

il. 

If the periodic orbit is stable or elliptic in nature 1I21I1 . then 
there is a finite area of phase-space surrounding the peri- 
odic orbit with very similar dynamics, i.e., an equivalent pre- 
dictable, periodic motion through phase-space, modulated by 
an approximately harmonic oscillation around the actual peri- 
odic orbit. The parameter regimes under which such dynamics 
occur can readily be determined by examining stroboscopic 
Poincare sections, produced by plotting the solutions of Eq. 
(I36> [modulo 2n in J'(J3)], for a large number of initial con- 
ditions and values of n, on the same axes. If, as in Fig. [2 
noticable island structures, centered on the points making up 
a stable periodic orbit, are apparent, then the action dynamics 
experienced by a significant fraction of the phase-space are 
approximately described by Eq. ( l37l . 
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FIG. 1: (color online). Poincare sections determined by Eq. <36> 
(white dots), superimposed on Wigner functions |34] corresponding 
to single wavepackets of the form given in Eq. Q (color density 
plots). Means and variances determined by Eq. <60> . for fc = 2 and 
(a) 6 = -0.2, (b) e = 0.2. Units are dimensionless. 



is, unlike J', explicitly dependent on the value of the quasi- 
momentum j6. 

If the initial phase-space distribution covers at least one 
phase-space cell [Itt x 2n in 6(J3) x J^(J3) space], then each of 
the islands surrounding a given stable periodic orbit are pop- 
ulated. We can then describe an averaged evolution of I„(J3). 
Furthermore, if the initial phase-space distribution is centered 
around I(j3) = and has a finite width in I(J3), with increas- 
ing n, the value of I(J3) can be considered progressively more 
negligible. 

Under these conditions, which are generally applicable to 
the atom-optics experiments carried out at Oxford 1 1, 2, 3, 4 
Hil^Qj, the momentum evolution of those initial conditions 
located within the island structures centered on a particular 
(p, j) periodic orbit is given approximately by 



I„(j3) - 2nn 



+ sgn(e)Q 



(39) 



As in this approximate formulation all jS-dependence of the 
evolution has been removed, it is a straightforward matter to 
produce an equivalent expression for the mean change in mo- 
mentum of cold atoms (in a freely falling frame) 



3. Prediction of quantum accelerator modes in atom optics 

With Eqs. j33> and ( I35> . Eq. (j^) can be transformed into 
an equivalent expression for I„p(J3). 



J„pO0) = I(/3) + 2nn[\ + sgn(e)pQ]. 



(38) 



As each of the /? subspaces is now in an equivalent, freely 
falling with gravity, y6-independent frame, the value of 



Pn - 27rn 



- + sgn(e)0 
P 



no 



(40) 



It is enhancements of the population of the observed mo- 
mentum distributions of cold cesium atoms, near these par- 
ticular values, that constitute quantum accelerator modes, as 
seen in the Oxford atom-optical experiments Ql |^ |3l |^ IH |^ 
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IV. SECOND-ORDER CUMULANT ANALYSIS OF 
QUANTUM ACCELERATOR MODES 

A. Derivation of the second-order cumulant map 

I. Overview 

The quantities we consider are the mean angle 9{j3) - 
{0}(J3), the mean action J(J3) - {J){J3), the angle vari- 
ance cr2(j6) = <0^>(y6) - Gifif, the action vaiiance S^(IS) = 
{u)(fi) - JTipy', and the symmetrized action-angle covari- 
ance T(y8) = {6J + je)(J3)/2 - 0{J5)J(J3), corresponding to 
the general quantities yu^, fi^, crl, crl, and crl., as described in 



Section lirXI respectively. 



5. Symmetrized action-angle covariance 

We first examine the covariance dynamics. FromEqs. j34at 
and (14 1> . we deduce the mapping equation for the covariance 
to be 



■r„+iOS) =^<[0„ + sgn(e)J-„+i]J-„+i>(/3) 

+ ^<J'«+i[0« + sgn(e)J-„+i])Oe) 
- mp) + sgn(e)J-„+i()8)]J'„+i()8). 

Substituting in Eqs. ( I34bt and (I42> . we get 



(46) 



2. Mean angle 

Taking the expectation value of Eq. ( l34aL we trivially de- 
duce 



0«+iOS) = 0«OS) + sgn(e)J-„+iOS). 



Mean action 



(41) 



Using Eq. JCl 1> to determine {smd„}(J3) - 
g-o";;(/3)/2 gjjj 6i„(;6) (derivations for this and a number of 
other useful identities are in Appendix |Cj, it is straightfor- 
ward to deduce from the expectation value of Eq. (I34b> the 
corresponding equation for the action variable: 

Jn^m ^Jn(J3)- le-^'-^^l^ sin 6„(J5) - sgn(e)27rQ, (42) 
which couples to the angle variance cr^ifi). 

4. Angle variance 

From Eqs. (|33a) and J41> . the corresponding mapping equa- 
tion for cr^(/6) is given by 

- 6„(J3f - 2sgn(e)6„(J3)J„^i(P) - Jn+xifif 
Substituting in Eqs. (|34b} and J42> . we get 

+ sgn(e)<§„ J-„ + J A - 2ke„ sin §„)(/3) 
- 20„(yS)[sgn(e)J-„(/3) - 'ke-'^'-^"^ sin0„C0)]. 



(43) 



Finally, making use of Eq. ( IC15> to determine (On sin 6n)(J^) ■ 
e-<^l(P)l-ie„(J3) sme„(J3) = crl(J3)e-'^'^^'^ cos e„(fi), we deduce 



2sgn(e)[r„(/3)-kcrf,(/3)e 



-o-m/2 



cos0„(J3)l 



(45) 



which couples to the action variance S^(J3), and to the sym- 
metrized action-angle covariance T(J3). 



rn^m =^{0„[J„ -ksm0„] + [J„ -ksme„]0„m 
+ sgn{e)Sl,(/3). 



(47) 



We once again use Eq. ( IC15> to evaluate {0„ sm6„)(J3), finally 
arriving at 



r„^m =T„(/3) - ~kcrl(J3)e-''"^''>''- cos e„(J3) 
+ sgn(e)5,ii(y8). 



(48) 



6. Action variance 

From Eqs. (I34b> and ( I42> . we deduce the mapping equation 
for the action variance to be 



Sl^i(/3) ={Jl + P sin^ 9 - ~k(J„ sin + sin 0J„))(/3) 
-J'«()S)'-feV""*^'sin2 0„(^) 
+ 2lcJ„(fil)e-'^'^^'^ sin 0„(I3). 



(49) 



We make use of Eq. JC13> to determine (sin^ 0(fi)} — 
e-'^' sin^ 0„P = {[1 - e"'^"][l + e"'"- cos 20(^8)] )/2, and 
Eq. (En|i to determine {J„(fil) sin 0(fil) + sin 9{J3) J ..(fi)) - 
J„(/3)e-'^"'^'>'^sin0„(P) = e-'^'^'J^'>'^r„(fi)cos0„{/3). Substitut- 
ing these results in, we finally produce 



^^^^ Sl^i(J3) ^Slifi) - 2ke-''"^^'^r„(fi)cos0„(/3) 



P 2 2 (50) 

+ 2-[l - e-'""<^^][l -H e-'^"'^>cos20„(je)]. 



7. Coupled mapping equations 

We thus conclude with five independent coupled mapping 
equations, approximately describing the 0(fi), J^ifi) kick-to- 
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kick operator dynamics in terms of means and variances: 

0„+ 1 (yS) =0„ + sgn(e) J-„+ iQ3), (51a) 
J',,+1 (yS) ^Jn(P) - h-'^"'^^'^ sin 0„(J3) - sgn(e)27rQ, (5 lb) 

+ 2sgn(e)[T„(/3) - fcfr^Oe)^--"*^'/^ ^^^^^^^ 



T„+i06) =Y„08) - ~kcrl(fi)e-''''^'>'^ cos 0„(;8) 

+ sgn(e)52^i08), 
5^^106) ^Slifi) - ih-'^'-'^^'^r „(J3) cos 6„(J3) 



(51c) 



(5 Id) 



P . 2 (51e) 

+ _[i +e-^„(/s) cos20„(;8)]. 



lated quantum accelerator modes, seemingly well explained 
in the e-classical picture of Eq. j36l l, at relatively large values 
of e. Note also that although the e-classical phase-space struc- 
ture determined by Eq. i36\ is clearly dependent on the value 
of k (or e^°~^^k), the prediction of the momentum evolution of 
those atoms occupying quantum accelerator modes [Eq. j40H 
is not. The area of the relevant island structure dictates the 
proportion of atoms to be accelerated, but not their averaged 
momentum evolution. This is further confirmation of the ro- 
bustness of quantum accelerator modes as an experimentally 
observable effect. 



2. Example: (p, j) = (1,0) quantum accelerator modes 



We note that setting all the second-order cumulants, cr ifi). 



S (J3), and T(J3), to zero causes Eqs. ( I51a> and ( I51b> to re 



vert to the effective classical mapping described by Eqs. ( I36a> 
and (I36b> . In the context of cumulant hierarchies, Eq. (I36> 
can be logically ordered as a first-order cumulant expansion 
of the dynamics of the underlying operator-valued variables, 
as described by Eq. ( I34> . The term "e-classical" is seen to be 
appropriate, as it is only in the limit of the commutators (pro- 
portional to e) tending to zero, that all quantum fluctuations 
can be summarily neglected. The second-order cumulant map 
described by Eq. (15 U is thus a lowest-order description that is 
still able to account for some of these effects. 



B. Periodic orbits in tlie second-order cumulant map 

}. Validity of the second-order cumulant map 

When considering under what circumstances we expect the 
second-order cumulant map of Eq. ( 15 1> to be a reasonably 
complete description of the dynamics, we note that Gaussian 
wavepackets, which are described completely by their means 
and variances, are stable when undergoing simple harmonic 
oscillator dynamics, and that in the vicinity of elliptic periodic 
orbits, the local dynamics are approximately harmonic i2Tll . 
We thus expect Eq. ( 15 U to be most useful when describing the 
dynamics of quantum accelerator modes, which is of course 
the situation of most interest to us! 

We are thus most interested in what occurs when the initial 
condition is located on top of an elliptic periodic orbit, largely 
inside a stable island of the e-classical phase space described 
by Eq. ( I36t . If the initial condition is located in the chaotic 
sea of the e-classical phase space, Eq. ( 15 11 1 is unlikely to have 
much predictive power, due to the more complicated underly- 
ing dynamics. 

We note that if the angle variance remains constant from 
kick to kick, i.e., cr^^j(^) = crliji) - a^, then Eqs. ( I51a> 
and ( I51b> again effectively revert to the e-classical map of Eq. 
( I36> . The only difference is that the kick strength is scaled by 
a Gaussian function of the angle variance, i.e., k is replaced 
by e^°"l^k. A stable solution to Eq. ( 15 1> would therefore ex- 
plain the experimentally observed presence of highly popu- 



Notable among these are the originally discovered (p, j) = 
(1,0) quantum accelerator modes around T — T1/2 111, cor- 
responding to a periodic orbit of order 1 and jumping index 
0, i.e., a fixed point, clearly observable at values of |e| up 
to ~ 0.8 |4]. We extend the usual procedure of determin- 
ing solutions of this kind in the e-classical map, where one 
sets 0„+i(j3) = = e, and J-„+i(j6) = J„(fil) = J, to im- 

posing these conditions, as well as cr\^y(fi) - (T^ifi) - <t^, 
T„+,(/3) = y„(fi) = T, and Sl^^iP) = Slip) = on Eq. 
This immediately produces 



sin§ = - sgn(e) — 



(52a) 
(52b) 



P ^- 2sgn(e)( Y - ka^e-^''^ cos 0), (52c) 
S 2 ^sgn(e)ka^e-'^''^ cos 0, (52d) 

2ke-^''^r COS0 ^—(1 - e-^'){ l + e-^\os20). (52e) 

We note that Eqs. ( I52a> and ( I52b> can be considered identical 
to the equivalent equations which would be produced when 
trying to determine fixed points in the e-classical map, with 
k replaced by e"°" ^^k. This follows directly from the similar 
result in the second-order cumulant mapping [Eqs. ( I51a> and 
JsTbll l. 

It is now desirable to reduce these coupled equations [apart 
from Eq. the fully reduced ( l52aH to closed form for the indi- 
vidual variables. We begin by equating Eqs. ( I52ct and ( I52dt . 
to determine that 



t = -a-V^'/^cosg. (53) 
Substituting this result into Eq. ( I52e> we get 



d-^e-^\l- sin^ 0) = - ( 1 - e-"^' ) [ 1 + e""^' ( 1 - 2 sin^ 0)] . (54) 



We now substitute Eq. ( I52b> into this, and order the result by 
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-2 

powers of e . We arrive at 



=e-^^' +2 5-2 



(55) 



a closed equation for the desired fixed-point angle variance 
(T^. With the aid of the quadratic formula, this can be factor- 
ized to 



= 



/27rQ\2 






(— )m 


(^2)2 + 

















/27rQ\2 






2 




(^2)2 + 






\-Tl-\ 







(56) 



Thus, one of the this product of two terms must equal zero for 
there to be a fixed point involving the variances as well as the 
mean values. That this can in fact never occur for finite tr^ is 
most easily shown graphically. 

The quantity (2;rQ/^)2 can be considered a freely varying 
parameter However, from Eq. ( I52bt . we can see that for there 
to be a hypothetical fixed point, then it is constrained by 



e"*" sin^g < 1, 



(57) 



as is clearly also true in the e-classical limit (e"°" — > 1). 

The quantity e""^' is compared with /-(o"^) - (2nD./k)^ - 
a-- - -\/(o^2)2 + [(27rO/fc)2 - 1]^ for various values of < 
(InQ/lcf < 1 by plotting both as functions of &^ in Fig.|3a). 
In Fig.|2jb) the same is done for f+{cr^) = (InQ/k)^ - cr^ + 
V(<t2)2 + [(27rQ/^)2 - 1]2. It can be cleai'ly seen that both 
f-id-^) and /+(o"2) only ever intersect with e^°" at cr^ - 0, 
which, as an extended solution cannot be considered reason- 
able in this context, is exactly the pointlike e-classical limit. 

The function e^°" is more generally approximated by 
/_(o"2) than f+{&^), where additionally {Inillk)^ must equal 
1, however. This can be seen by expanding f-{&^) as a 
McLaurin series: 





FIG. 2: Comparison of e 
(T^ - V(o^^)^ + [(2;r£l/^)2 ■ 



(solid line), with (a) /_(<t2) = {InSljkf 
and (b) f^ia'-) = {InSllkf - tr^ 



V(<t2)2 + [(27ra/fc)2 - 1]2, for values of (InSl/kf equal to (points), 
0.2 (circles), 0.4 (crosses), 0.6 (plusses), 0.8 (asterisks), and 1 
(squares). Note that /_ = e"""" only when a- = 0, and that /+ = e"'^^ 
only when cr^ = and (InCl/k)^ = 1. 



C. Constraining the second-order cumulant map to be 
uncertainty conserving 

1. Motivation 

Our result on the impossibility of fixed points in the second- 
order cumulant map [Eq. ( I51H . can be understood with the 
aid of some simple considerations. Recall that finding a 
fixed point of e(J3), J(J3), o-^(J3), T(y8), and S^iJS), is equiv- 
alent to finding a perfectly Gaussian kick-to-kick eigenstate 
of the system. Even if one has an initially perfect Gaussian 
wavepacket, fully described in terms of its means and vari- 
ances, it is not possible for it to remain Gaussian unless no 
terms in the Hamiltonian are of greater than quadratic order 
in the observables, or equivalently, that no terms in the corre- 
sponding Heisenberg equations of motion are of greater than 
linear order. An actual Gaussian eigenstate can only occur for 
the harmonic oscillator 

Although the local dynamics in the vicinity of an elliptic 
periodic orbit can be considered approximately harmonic, the 
global dynamics described by both the e-classical map [Eq. 
( I36H and the underlying Heisenberg map [Eq. ( I34H are in fact 
highly nonlinear in the observables. Already at second order, 
the cumulant dynamics accurately reflect this fact, and this 
will be true for dynamics occuring near any e-classical stable 
periodic orbit. 



^2n2 



Md-^) =1 - + 



2[{2nQ.rkf 



1] 



3[(2;rQ/;t)2 - 1]3 



(58) 



which is equal to e °" to first order in and to second order 
in cr^ if 2nQ.lk - 0. Thus, as shown in Fig. 13 a), the best 
matching of f-(&^) to e"°~ occurs near &^ - 0, and when the 
effect of gravity is relatively small. 



2. Enforcing uncertainty conservation 

Clearly to find solutions which are approximately stable 
fixed points in terms of the first- and second-order cumulants, 
we need a way of constraining the dynamics such that an ini- 
tial perfectly Gaussian wavepacket remains Gaussian after ev- 
ery iteration. 

Recalling the generaUzed relation for Gaussian wavepack- 
ets given in Eq. Q, we can state that for a Gaussian 
wavepacket (r2^;(;8)52^,(^)-T„+i(y8)2 = e2/4. Replacing Eq. 



11 



( 15 let with this constraint forcibly slaves the dynamics of the 
action variance to those of the position variance and the sym- 
metrized action-angle covariance. In the present context, this 
appears to be the simplest way to enforce genuinely Gaussian 
dynamics. 

We substitute Eqs. dsTcl and dsTdl into (rl^i(fi)Sl^^(J3) - 
T„+i(J3)^ = and produce, after some straightforward ma- 
nipulation: 



[r„(/3) - kcrl(J3)e 



COS 0„(/3)Y + e^/ 4 



(59) 



Equation ( I59> can in turn be substituted back into Eqs. ( 15 let 
and ( I51d> . and we end up with a closed set of four coupled 
mapping equations: 



Jn+m ^J,m - fee-""'^'/' sin e„(J3) - sgn(e)27rn, 

cr^^iOS) ^cr^OS) + 2sgn(e)[Y„08) - ~ke-''"^^'Wjfi) cos 6„m + 



[T„Og) - ^e-".^'^>/VgOg) cos 0„m^ + eV4 



T„+iOe) =T„Oe) - ~ke-''"^^'^crl(J3) COS 9„Q3) + sgn(e) 



[r„Q3) - fce-^'^>/Vg(/j)cosg„Og)]^ + 6^4 



(60a) 
(60b) 

(60c) 
(60d) 



where S^^^(J3), if desired, can be deduced from 
^l+i(fi)s£i(f^) - T„+i06)2 = e2/4. We have reduced to 
four independent equations in Eq. ( I60> from the five of Eq. 
(15 U . due to the fact that the evolution of the action variance 
S ^ifi) is no longer independent. 



Fixed points 



We once again impose the conditions 6„+i{J3) - 0„(p) - 6, 
Jn^i(j3) = J„(J3) = J, (rl^.ili) = (rlili) = d-\ and Y„+i06) = 



f„(J3) - Y, this time on Eq. ( I60> . in order to determine any 
fixed points. The condition S^^^^(j3) - S^ifi) = is thus 
fulfilled automatically. 

The result of imposing these conditions is given by: 



sin 6 ■■ 



sgn(e)- 



2t ^2ke-^'^^a^cose 

[t - ke-'^'l^&^ cos ~6f + e^lA 



(61a) 
(61b) 



-sgn(e)- 



(61c) 



ke 



(T COS 6 =sgn(e) 



[r -ke-^^^l^a^cos 0]^ + e^/4 



(61d) 



where we note that the conditions for the mean values [Eqs. 
(I61a> and ( l61bH are exactly the same as those produced in the 
case of the full second-order cumulant map [Eqs. ( 15 2 at and 
( l52bH . and thus the corresponding equations in the e-classical 
limit, with k replaced by e^°"^-k. 



When reducing the equations for the variances to closed 
form, we expect differences to occur in the equations for the 
variances, however, although adding Eqs. ( I61c> and ( I61d> to- 
gether yields Y = ke^°"^^(r^ cos 6/2, an identical result to Eq. 

Substituting this back into Eq. ( I61d> yields 



Y^ + eV4 



sgn(6)Y 



(62) 



which reveals that sgn(e)Y must be positive. 

Combining the squares of Eqs. iblh\ and ( I54L along with 
the positivity of sgn(e)Y, yields a second independent equa- 
tion dependent only on cr^ and Y: 



sgn(e)Y : 



Ikd-^e-'^'/^ 



(7rQo-2)2 



(63) 



Substituting this back into Eq. i62\ produces 



2cj2 



(7rQo-2)2 



(kiT^e-^'l^\ 



-(nna^f + — , 



(64) 



a closed equation of 5" , for which, depending on the values 
of k, 2nQ., and e, solutions do in fact exist. Such a Gaussian 
solution is thus dependent on e as well as k and 2jtQ.. 

This equation can therefore be used as a starting point 
to determine fixed points, for the effectively Gaussian map- 
ping of Eq. ( I60> . numerically. For the fixed points stud- 
ied in this Paper, we consider situations which correspond 




E 



FIG. 3: (color online). Number of iterations of Eq. <51> evolved by 
a Gaussian stable fixed point such that \J'\ < n. Black indicates ab- 
sence of 6-classical stable solutions 1 18], white absence of Gaussian 
stable solutions. Numbers label the contours where \J'\ < n for that 
number of iterations (the number of iterations is capped at 100). The 
solid line marks the average experimental laser intensity 0^ = 0.8;r, 
dashes demarcate its experimental range (0.3;r-1.27r) |;4||. Units are 
dimensionless. 



experimentally to freely varying the kicking periodicity and 
the laser intensity, with InQ. - gGT^ determined hy T - 
(e + 2n)Tii2l2n, g = 9.8 ms^^, and G = 27r/(447 nm). 
For illustrative purposes, Wigner representations [ W(9, J') - 
(27r|e|)-i / dTe-'^^'^'^i/j'{9-T/2)if,{0 + T/2) ] jH of such "sta- 
ble" Gaussian wavepackets, overlaid by Poincare sections of 
the e-classical phase space lUsll . are shown in Fig.^ We see 
that the Wigner functions closely match the shape of the stable 
island Ii34ll . 

It should be noted that this same general procedure for find- 
ing approximate fixed points can be applied, with slight mod- 
ifications, to general periodic orbits. A simple transformation 
of the action, i.e, defining X.n(J3) = J^niP) — 27rni/p, causes a 
periodic orbit, of order p and jumping index j in Oifi), J'(J3) 
space, to have jumping index in 0(J5), X.(J3) space. One can 
then search for solutions of the resulting map, iterated p times, 
in an analogous manner 



D. Exact wavepacket dynamics: Comparison 

We now come to the final point of our analysis, propagat- 
ing Gaussian fixed point solutions using the full second-order 
cumulant mapping, where all variances must be considered 
explicitly. 

In Fig. |3] we display the time for which the center of 
mass momentum remains inside its initial phase space cell 
{\J'„(J3)\ < n), using this as a rule-of-thumb measure of rela- 
tive longevity (for a genuine fixed point this would be forever). 
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FIG. 4: (color online). As in Fig.|3| but for exact wavepacket evo- 
lutions, propagating initial conditions from Eq. <67t with the time- 
evolution operator of Eq. {8} (/3 = 0). Units are dimensionless. 

We see that there is a sizable region where there there are no 
stable Gaussian fixed points, in addition to the region where 
there are no e-classical solutions, and that quantum acceler- 
ator modes for e < are generally more long-lived. These 
observations are broadly born out by experiment (H]. 

We have also computed for how long \J'„\ < n when inte- 
grating the exact evolution described by Eq. ©IS]. Due to 
the fact that we restrict ourselves to a single /3 subspace, the 
wavefunction we propagate must be a Bloch state, of the form 

{Gz = 2;r/ + 6\t//) = e''^'-^'''^^'>>J/(e), (65) 

which has a discrete momentum spectrum. As we wish to 
carry out the numerics in the momentum basis, we are there- 
fore obliged for practical reasons to work in the I(J3) basis 
rather than that of J'. The reason is that the spectrum of 
I - \e\k is fixed, wheras that of I moves with time [see Eq. 
j33lll. As the discreteness scale is |e|, this issue vanishes in the 
e-classical limit, but it is unavoidable for any finite value of e 
in the fully quantum-mechanical dynamics. 

Determining the appropriate initial value for I(J3) from a 
predetermined value of is, with the aid of the expectation 
value of Eq. J33l l. a straightforward matter. The values of 
and the variances 5"^, T, and S^, are unaffected by this trans- 
formation, and so we can set the initial i^{0) to be 

if/(6,/3) oc[27ro-2]-'/4 

/ [l-/2t/|6|][g-0]^ iI(MO-0] \ 

(66) 

The computationally more convenient I representation of 
the initial state determined by Eq. (16 1> is then the discrete 
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Fourier transform of Eq. 
n|e|), where 



l66> . given by [i/r) oc 2^-coC«|J - following operator- valued kick-to-kick map: 



c„ = 



<r2 



1/4 



27r(e2/4 + T2)_ 

1 - /2t/|e| 



X exp - 



fe[n\e\-im 



(67) 



This must in principle be normalized numerically, although 
for a Gaussian state which is well localized in 6, i.e., one 
where Vcr^ <K In, the effect of this will be negUgible. 

Figurel^shows the results of these integrations. We see that 
Fig-Ereproduces its qualitative features quite well, especially 
for smaller values of e and k. More surprising is the repli- 
cation of a saddle-point feature at around {e - -1.5,^ = 2), 
indicating a resurgence of stability for large e that is clearly 
not an artefact of our approximations. 



V. CONCLUSIONS 

In conclusion, we have developed a general method for 
using second order cumulants to study semiclassical-like 
dynamics near stable periodic orbits in phase space. We 
have successfully applied this method to quantum accelerator 
mode dynamics, which operate in an unusual e-semiclassical 
regime, thus gaining insight into the longevity of quantum ac- 
celerator modes in different parameter regimes. We have ex- 
plictly determined the second-order cumulant mapping for the 
dynamics of the quantum 5-kicked accelerator, taking place 
near integer multiples of the the half-Talbot time. In the pro- 
cess, we have shown the e-classical dynamics derived by Fish- 
man, Guarneri, and Rebuzzini, to be the first order approxi- 
mation within the relevant cumulant hierarchy, we have ex- 
plained why this description remains effective even for non- 
negligible values of e, when describing quantum accelerator 
mode dynamics, and we have shown why, as quantum acceler- 
ator modes are traced back to stable periodic orbits in the rel- 
evant phase space, our methodology is particularly well suited 
for analyzing quantum acclerator mode dynamics. 



T gGT^ 

Zn+l ^Zn + —Pn+\ (Ala) 

m 2 
Pn+i =Pn - Gh4>d sin(Gz„) + mgT. (Alb) 

In order to reduce the number of free parameters to a mini- 
mum, we rescale the dynamical variables to be dimensionless, 
such that y = Gz and p = GTp/m. This produces, from Eq. 

(ED EH, 



Pn+1 -Pn - K?,mXn - 27rQ, 



(A2a) 
(A2b) 



where k - G^Tfupd/m is the usual dimensionless stochasticity 
parameter associated with the (gravity-free) classical 5-kicked 
rotor |21], and, as elsewhere in this paper, 2nQ. - gG^T is a 
dimensionless parameter quantifying the effect of gravity. 

The conventional semiclassical limit of these equations can 
be achieved by replacing the operator-valued dynamical vari- 
ables, X p, with their expectation values, x - ix) and 
p = ip). As explained in Section IIV Al this must coincide 
with the vanishing of the commutator 



-i\X,p\ =k = 



HG^T 2nT 



1/2 



(A3) 



where ri/2 — Inm/tiG is the half-Talbot time. 

The scaled Planck constant k is that conventionally used in 
studies of 5-kicked rotor/particle dynamics, particularly when 
investigating the semiclassical limit, defined as 
We recall that the smallness parameter e upon which the con- 
cept of e classics and semiclassics in this Paper is based, is 
defined within 



r = (2;r^+e)^ = r,/2 



(A4) 



where £ eZ. Substituting Eq. (IA4> into Eq. ( IA3t . we therefore 
get 



k^27i{ + e. 



(A5) 
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and we see that taking the e-classical limit is equivalent to 
the conventional rescaled Planck constant k being equal to an 
integer multiple of 2n. 



APPENDIX B: BLOCH PARAMETRIZATION OF THE 
DYNAMICAL VARIABLES 

1. Overview 



APPENDIX A: CONNECTION BETWEEN e AND THE 
CONVENTIONAL SEMICLASSICAL LIMIT 

From Eq. (|5}, the Heisenberg equations of motion for z and 
p can be determined directly, and integrated to produce the 



We will now explore some relevant details of our chosen 
parametrization of the position and momentum variables, for- 
mulating inner products, operators, and commutation rela- 
tions in terms of the angle 9, the quasimomentum yS, the dis- 
crete position /, and the discrete momentum k. 
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2. Inner products 



3. Operators 



We define 



1 



V27r 



and, bearing in mind the limited ranges of e [0, 2ti) and 
yS e [0, 1), note the following identities: 



-ik{6-ff) 



--2n6{0 - ff). 



k=-co 



J dOe" 



■e(k-k') 



-2n6kk', 



l=—oo 



--dw. 



(B2) 
(B3) 
(B4) 
(B5) 



Inserting the position representation of the identity operator 
into the inner product {Gz - 2nl + 6\Gz - 2nl' +6'} we find, 
using Eq. (IB1> . that 



(Gz = 2nl + e\Gz = 27tI 



^ y 

2n 



%0-0') 



(B6) 



k=—oo 



Substituting Eq. (IB2> into Eq. JB6> . it simplifies to 



{Gz^27il + 0\Gz^2nl' + &) ^ J d^e''^^''^'-''^6(0- 6'), (B7) 



which, with the aid of Eq. ( IB3> . results in the following gen- 
eral form for the inner product of two rescaled position eigen- 
states, in terms of the angle 6 and the discrete position variable 



/: 



{Gz = 2kI + 6\Gz = 27:1' + 6') = 5u'6{e - ff). 



(B8) 



Similarly, using Eqs. ( IB4t and ( IB5t . one can straightfor- 
wardly determine a corresponding expression for the inner 
product of two rescaled momentum eigenstates, in terms of 
the quasimomentum yS and the discrete momentum k: 

{(tiGY'p ^k+/3\(nG)-'p = k' +15') = 5kk'5(Ji-l3'). (B9) 



Having determined all the relevant inner products in terms 
of 0, 1, j3, and k, we set 



/oo 
dO^jlGz^ 2kI + e)l{Gz = 2nl + 0\, (BIO) 
/=-co 

/ oo 
de \Gz = 27rZ + e)9{Gz = 27rZ -H 6I|, (B 1 1) 

/=-co 

/ oo 
£/yS ^ \(nGr^p ^k+ l3)k{{nGy^p ^k+pi (B12) 
jt=-co 

/ oo 
^/yS 2] K^G)-^ = ^ +P)mGr'p = ^ (B13) 
*:=-co 

These operators are defined such that l\Gz - 2nl + 9) - 
l\Gz = 27il + 0), 6\Gz = 27TI + 0) = e\Gz = 27rZH-6'), ^K^iG)-'/? = 
k+P) = k\ihG)-^p = /t -H /?), and MhGy^p = yt + /?) = 
/5\{hGy^p - k + P), as can readily be confirmed with the aid 
of Eqs. (IB8> and ( IB9> . In addition, clearly 

+ 27tI^Gz, (B14) 

P + k^^. (B15) 
«G 



4. Commutators 

We now determine the commutators of the various 
parametrized position and momentum dynamical variables. 

We begin by considering the angle 0, and the quasimomen- 
tum p. Using Eqs. JB11> . jB13> . and (IB1> . we readily deter- 
mine that 



pe=-^ ff dpdepoe 

V27r J J 



-iB/3 



X \{nGr^p^k+P)e-'^'"'^^''^^\Gz^2nl + 9\. 



k.l=-a 



(B16) 



Sandwiching this expression between identity operators, in 
the position and momentum representations, respectively, pro- 
duces, with the aid of Eqs. (jED, and (jUJ, 

pe=-^ ff dpdopoe'^^ 

V27r J J 

OO 

X \(^G)-^p = k+P)e''"''e'^'''''^{Gz^27Tl + 9\ 

k' J'-—oo 

=(y60)"'" 

^ep. 

(B17) 

It therefore follows that the quasimomentum and angle oper- 
ators commute, i.e.. 



{e,p] = 0. 



(B18) 
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Considering now the discrete momentum k and the angle 9, 
we determine, from Eqs. JBH . JBl 1> . and JB12> . that 



kO-- 



k6e 



-iOk 



d/3 Yj K^GT^p = /fc + ^ 27r/ + 6\. 

l=-oo 

(B19) 

In analogous fashion to what was carried out for the pO prod- 
uct, we sandwich this expression between identity operators 
in the position and momentum representations, respectively, 
producing with the aid of Eq. (IB1> : 



1 rr °° 
ko I dp'de' y \(fiGy^p = r +/?' 

/oo 
d60e'''' J] /te'' 



i\P' (ff +1kI' )] 



We now use the identity 



'k(e'-e) 



(B20) 



dem^6{e-e')^^jie') 



dff 86' 
along with Eq. JB2l i to determine that 



(B21) 



J 00' 



-0) 

k=-cc 

= -;27r I d09e" "^6(0-6') 
80' 

=2n{ffk' - Oe"^"' 



(B22) 

and insert this expression back into Eq. jB20l l. to reveal 

k0 = (B)^ -i^0k- i, (B23) 

which result is more conveniently phrased within the commu- 
tator 

[0, k] = i. (B24) 

In a very similar manner, using Eq. JB4> and inserting a 
partial derivative with respect to the quasimomentum /3, it can 
be shown that Inl^ = (2nip)'^ + i, and therefore that 



[/,jS] = 



27r 



(B25) 



Referring to Eqs. JB14> and ( IB15l l, we can now deduce 
from [z,p]/n = [0,P] + [0,k] + 27r([[,yS] + [l,k]) = ' that 



[k,l] = 



2k' 



(B26) 



Trivially, [0, /] - [k,0] - 0, and so, to summarize: 

[0, k] = i, (B27) 

(B28) 

[0,;e] = [0,/]-[^,;6]-O. (B29) 



[k,l]^[l,/3]^—. 
In 



APPENDIX C: SECOND-ORDER CUMULANTS 
1. Counting statistics 

We choose to approximately express a given order 2n mo- 
ment {qiqi ■ ■ ■ qin) in terms of first and second order cumu- 
lants only, i.e., means and variances, setting third- and higher- 
order cumulants to zero [see Eq. ©]. The moment can be 
expressed as a sum of n terms, each consisting of a product of 
k variances and 2(n - k) means. 

Taking one such term, when considering how many ways 
there are to produce a product of k variances and 2{n — k) 
means, this is equivalent to determining how many ways there 
are to partition the terms into two subsets, one with 2k 
elements which pairwise form variances, and the other with 
2(n - k) elements individually forming means. 

Our starting point is to consider all possible permutations 
(i.e., orderings) of the terms. For n = 3 one possible per- 
mutation is 



q^qiqeqiqAqi- 



(CD 



The first 2k terms in the permutation are partitioned into 
unique neighbouring pairs, and all subsequent terms are are 
partitioned into individual units, e.g. for « = 4, A: = 2 and this 
permutation 



«§3§1»«§6§2»«§4»«§5»- 



(C2) 



We thus have k pairs, identified with variances, and 2(n - k) 
individual terms, identified with means. There are (2n)! such 
permutations. 

As there is only one "correct" ordering, set by the ordering 
of the operators in the original moment, we must correct for 
overcounting. In particular the "correct" ordering of the op- 
erators making up the variances is in ascending order of the 
value of the subscript. The ordering of the individual means 
and variances with respect to each other is obviously unim- 
portant. The complete set of equivalent orderings is thus: 

{«§1§3»«§2§6»«§4»«§5», 
«§1§3»«§2§6»«§5»«§4», 
«§2§6))«§1§3»«§4»«§5», 
«§2§6))«§1§3»«§5»«§4», 
«§3§1»«§2§6»«§4))«§5», 
«§3§1»«§2§6»«§5»«§4», 
«§2§6»«§3§1 »«§4»«§5», 
«§2§6»«§3§1»«95»«§4», 
«§1§3))«§6§2»«^4»«§5», 
«§1§3»«§6§2»«§5»«§4», 
«§6§2»«§1§3»«§4»«§5», 
«§6§2))«§1§3»«§5»«§4», 
«^3§1»«§6§2»«94»«§5», 
«§3§1»«^6§2))«^5»«§4», 
«§6§6))«§3§l))«§4))«§5», 

«§6§2»«§3§1»«§5»«§4))) ■ 



(C3) 
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In general there are k\ (all the possible orderings of the vari- 
ances) multiplied by 2* (all the possible operator orderings 
within the variances to the power of the number of variances) 
multiplied by [2(n - k)]\ (all the possible orderings of the 
means) equivalent possibihties. Only one should be counted. 



Drawing on Eq. (IC4> it is now straightforward to approxi- 
mately expand the expectation value of the cosine of an oper- 
ator in terms of means and variances: 



2. General expansions of moments in means and variances 



If the qi are identical, then, designating the mean {{§)) = 
(q) - fi and variance {{(f)) - {q') - {q}^ - cr^, we can 
determine an approximate expression for the moment {q^") 
purely in terms of fj.^ and cr^: 



l^[2(«- 



(2«)! 



[2(«-fe)]!A:!2* 



(C4) 



For a general order 2n -H 1 moment, the corresponding approx- 
imate expression is given by: 



2„,1 ^ y (2n+l)\ ^ 2(«-*)+l( 2)*^ (C5) 



If we consider the expectation value of an averaged sym- 
metrized sum of products of a separately considered opera- 
tor q[ with a product of identical operators qj, i.e., {qiq^" + 
qfqi}/2, then Eqs. dJ and O can be combined to produce 
a similar approximate expression: 



(2«) 



k=0 



P2 



(2n - 1)! 



+ 2ncr?T 

12 [2(n-k)- lV.kl2^ 



(C6) 



x/.f-*^-'(cr2)\ 

where yUi - (qi) and fj.2 = (§2) are means, a2 - {q^ — {qif' 
is a variance, and crn - {qiqi + §2§i)/2 - (qOiqi) is a sym- 
metrized covariance. Similarly, for an equivalent symmetrized 
moment containing odd powers of q2, 



(2n-l- 1)! 



[2(n-k)+ l]!;t!2* 



+ (2n+l)o-]^ y !- 



(2«)! 



[2in - k)]\kl2'' 



(C7) 



{cosq)^y^-^{q'") 



(2«)! 
(-1)" 



l^(2n)\{-^A2{n-K)\\k\2^^ 



«=0 ^ ' k={) 



V' V (1)" 2(«-*:)/ 2 /ON* 

h^[2(n-kW.kl^ ^"""-^ 



k=Q n=k 



(C8) 



(-cr^/2f (-1)' 2, 



/t=0 /=o 

=0 

cr=/2 



z 



—e ' cos //. 

One can similarly draw on Eq. dCSI l to carry out an equiva- 
lent expansion for the expectation value of the sine of an op- 
erator: 

(-1)" 



<sing)^y ^ ' <g^"+') 



«=o 

00 00 



ZjZj \2{n-k)+ \\\k 



2(„-.Hi(^/2)^ 



/t=0 «=*■ 



(C9) 



-fH -H* (2Z H- 



(2ZH- 1)!A;! 



yt=0 /=0 

-e^""^"^ sin/v. 

Finally, by making use of Eq. (IC7l i. one can straightfor- 
wardly determine an approximate expansion of the expecta- 
tion value of an averaged symmetrized sum of products of the 
the sine of one operator multiplied by a second operator, in 
terms of means and variances: 
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-<§isin§2 + (sin§2)§i) (2« ^ i ^' +^2"^'^'^ 



«=o 



(2n + l)!2 
(-1)" 



V-i/ yj (2»+l)! 
~ Zj (2n + 1)! Zj 1^' [2(« - A:) + 1]!A;!2*^ 



2(«-/:)+1 . 2^t 2 (2« + 1)! 2(„-k). 2.k 

^2 <^o-2^ +^i2ro^„_ 7Ani7.nt^2 "^"-jj 



00 CO 



4=0 *:=« 
" (-1/ 



(-1)" 



[2(n-;t) + l]!;t!^2 2/ ^^[2(n-k)V.k\ 



•[2(«-A:)]!fe!2*' 



= yizii(^2/2)^Liy^ 

Z. ^ 2/ ; Y^Z.(2U 



( 1)' ..2/+I , _2 ^J2',.2/ 

(20! 



1)! 



;={) 



A'2 



(CIO) 



With these identities we now have everything necessary to 
determine the second-order cumulant dynamics quantum ac- 
celerator modes in the ^-kicked accelerator. 



3. Formulation of expansions relevant to quantum accelerator 
mode dynamics 

The dynamical variables relevant to (5-kicked accelerator 
and quantum accelerator mode dynamics are 6 and J'. All 
expectation values are evaluated within a reduced subspace 
particular to a specific value of the quasimomentum /3, only, 
as detailed in Section lTll B 61 

Bearing this in mind, we determine from Eq. iC9i that 



<sin0)(yS) = e^^'i^'/^sineOS). 
Using Eq. (IC8> . we determine that 

<cos0)(yS) = e-'^''^'''^cos0(fi), 

and that 

(sin2 0)(yS)=i[l-<cos(20))(j8)] 

= l{l_e-2-Wcos[20(^)]). 



(Cll) 
(C12) 

(CI 3) 



Using Eq. ( IC10> . we straightforwardly determine that 

sin e + (sin 0)J)(j3) ^e-''^P'>I^J(j3) sin 6(li) 
2 (C14) 



contained within any such subspace is periodic, multiplied by 
a quasimomentum and position dependent phase e"'^*^^, and 
can be equivalently represented by a rotor wavefunction. We 
will now discuss the connection in some detail, in particular 
so as to more explicitly link the operator oriented formalism 
used in this paper with the y6-rotor picture used originally by 
Fishman, Guarneri, and Rebuzzini 1118(1 . 



2. Density operator 

We consider an assumed well defined and unit trace den- 
sity operator, which, in the momentum representation, has the 
general form 



k,k'=-oo 

X \inGr^p /3}{(nGr^p ^k' + p'y 



(Dl) 



We select out that part of the density operator specific to a par- 
ticular y6 subspace by sandwiching it between yS-specific pro- 
jection operators Pip), as defined in Eq. J28> . The resulting 
expression is given by 



+ e-''"<fi^'^r(/3)cose(J3), 'P(J3)p'P(J3)= d/3"d/3"' \ dk>k"d(/3 - /3") 



and, again by using Eq. (IC10> . where we set q\ - q2, we 
determine that 

0smO}(P) = e-'^''J^'''\e(J3)sm0(j3) + o-^(j3) cos 0(j3)]. (C15) 



APPENDIX D: ANALOGY: /3-ROTORS 
1. Overview 

Subspaces of different /3 (in an accelerating frame) are de- 
coupled in 5-kicked accelerator dynamics. A wavefunction 



k,k' ,k" ,k"'=-co 

xdk"k'"(/3",/3"')Sk"'k'6(JS"' -fS) 

X mcr^p l3){(nG)-^p ^k' +I3\ (D2) 

00 

= Yj dw 03,13) 

k,k'=—oc 

X mcy^p ^k+ pXihCr^p ^k' + 13\. 

A general matrix element of this projected out density opera- 
tor, in the position representation, is then given by the follow- 
ing inner product: 
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^ {Gz = 2nl + 6\(hGr^p = k+ P)dkk'(J3,j3){inG)-^p = k' + /3\Gz = 2nl + 6}^ g(e, 6' ,fi)e''^^'^-'''^e' 

I 



'2n(l-l')i 



(D3) 



k.k'=- 



where 



i(ke-k'e') 



k,k' 



In 



(D4) 



Just as for plane waves, states resulting from the projection 
P(J3)pP(J3) are extended in position space, due to the fact that 
they are, in general, incoherent superpositions of Bloch states 
for one particular value of the quasimomentum /3. 

Bloch states are not normalizable (i.e., they lie outside the 
Hilbert space) when the inner product is defined in the usual 
way; the trace of the expression Eq. ( ID3> , defined as summing 
and integrating over all terms where 6-6' and I - /', can be 
seen to be similarly problematical. 

We note, however, that the normalizing factor N(J3) defined 
in Eq. ( I27> is given by integrating over the diagonal elements 
of g{d, 6') alone, i.e.. 



r d6g(6,6,/3)^ dkk'(fi,li) f d6 



,i(k-k')9 



2n 



= 2 du(l3,/3) = 7V08) 



(D5) 



k=-oo 



This takes advantage of the fact that integrating the diagonal 
matrix elements of Eq. (ID3> over any length In interval will 



give the same answer This is in turn due to the periodic nature 
of the state, apart from quasimomentum dependent phases, 
which are in any case absent on the diagonal. 

Effectively, if we are restricted to a single /5 subspace, we 
can define an inner product as the integral over 6 from a to 
a + 2n, at a particular value of /, which can be chosen arbitrar- 
ily. This provides a well-defined norm for Bloch states, which 
have been mapped from the inifinite, continuous position ba- 
sis appropriate for a free particle, to the periodic position basis 
appropriate for a rotor 



3. Discrete momentum operator 



We consider the action of the discrete momentum operator 
k restricted to a particular/? subspace, i.e., when acting on a 
projected density operator as defined in Eq. (ID2> . 

In the position representation, we end up with the familiar 
differential operator form appropriate for dynamics occuring 
on a circle: 



{Gz = 2nl + 6\kP{fi)pPmGz = 2nl' + 6') =e'-2'^('-n/3g/(e-«')/J J_ V kdkk,(/3,j3) 



2n 



k,k'-—oo 



(D6) 



-i—g(6,6',B) 
86^ ^ 



Note that the differential operator -id 1 86 needs to be taken 
inside the quasimomentum dependent phases, and thus it acts 
directly on g{d, 6'). 



4. /^-conditional expectation values 



a. Discrete momentum expectation value 



k multiplied by the density operator p, i.e.. 



1 °° 

<^>(J^) = 777^ Z ^^^GT'p = k+p\kp\{nG)-'p ^k+p). 
N(J3) ^ 



k=—oo 



(D7) 

Inserting the general density operator defined in Eq. ( IDll i pro- 
duces 



The (normalized) expectation value of the discrete momen- 
tum operator k conditioned to a single value of /3 is given by 
taking the normalized partial trace, as defined in Eq. ( I27> . of 



J oo 



(D8) 
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With the aid of Eq. dD4t . some fairly elementary manipula- 
tions prove this expression to be fully equivalent to 



(km 



If °° 



, e 



i(k-k')e' 



da,(J3,/3)k- 



k,k'=-o( 



In 



k,k'=-co 



dkk'(P,/S) — 
2n 



j(ke-k'e') 



de 



= I d6\-i—md\/i) 



where [taking Eq. {D5\ into consideration] 



Nip) Jdeg(0,e,l3) 



(D9) 



(DIO) 



is the angle-dependent component of the y6-projected density 
operator of Eq. ( ID2t . normalized to the proportion of popula- 
tion present in that y6-subspace. 



b. Angle expectation value 



We see that, following some simple manipulations, we have 
arrived at an equivalent expression to Eq. ( ID9> for the angle 
variable, written in terms of an integral involving q{9, 6',/3), as 
defined in Eq. JD10> . only. 



General expectation values 



In a similar fashion to the derivations given in Eqs. ( ID9l l 
and (ID12t . it is now straightforward to deduce an equivalent 
general expression for /^-conditional expectation values: 



The corresponding expectation value for 9 conditioned to a 
single value of (3 is, similarly, 

^ oo 



k=—o 



(DID 

by definition. Once more we substitute in the density operator 
of Eq. iPH . as well as inserting the position representation 
form of the identity operator. 
This produces 



1 rr °° 

^^^^^ ^7m JJ '^^'''^ ? (^dk^kiP'^P) 



: k + P\GZ 



2nl + 0) 

X {Gz = 2nl + 0\{hGr'p ^k' + /3'} 

i(k-k')6 



(D12) 



\ c 



deer 



6' =6 



(D13) 



where alternative orderings can be determined with the aid of 
the commutation relations. 



We thus see that for all possible expectation values of inter- 
est, g{9, 9' ,13) is a density operator containing all the informa- 
tion we need. The dynamics of each y6-subspace can thus be 
mapped onto the dynamics of separate rotors, termed yS-rotors 
by Fishman, Guarneri, and Rebuzzini Lldl . 
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